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BK Ca 2+ -activated K + currents exhibit diverse properties across tissues. The functional variation in voltage- and Ca 2+ - 
dependent gating underlying this diversity arises from multiple mechanisms, including alternate splicing of 
Kcnmal, the gene encoding the pore-forming (a) subunit of the BK channel, phosphorylation of ot subunits, and 
inclusion of (3 subunits in channel complexes. To address the interplay of these mechanisms in the regulation of 
BK currents, two native splice variants, BKo and BK SRK r, were cloned from a tissue that exhibits dynamic daily ex- 
pression of BK channel, the central circadian pacemaker in the suprachiasmatic nucleus (SCN) of mouse hypo- 
thalamus. The BKo and BK SRK r variants differed by the inclusion of a four-amino acid alternate exon at splice site 
1 (SRKR) , which showed increased expression during the day. The functional properties of the variants were inves- 
tigated in HEK293 cells using standard voltage-clamp protocols. Compared with BKo, BK S rkr currents had a signifi- 
cantly right-shifted conductance-voltage (G-V) relationship across a range of Ca 2+ concentrations, slower activation, 
and faster deactivation. These effects were dependent on the phosphorylation state of S642, a serine residue within 
the constitutive exon immediately preceding the SRKR insert. Coexpression of the neuronal [34 subunit slowed 
gating kinetics and shifted the G-V relationship in a Ca 2+ -dependent manner, enhancing the functional differences 
between the variants. Next, using native action potential (AP) command waveforms recorded from SCN to elicit 
BK currents, we found that these splice variant differences persist under dynamic activation conditions in physio- 
logical ionic concentrations. AP-induced currents from BK SRK r channels were significantly reduced compared with 
BKo, an effect that was maintained with coexpression of the (34 subunit but abolished by the mutation of S642. 
These results demonstrate a novel mechanism for reducing BK current activation under reconstituted physiologi- 
cal conditions, and further suggest that S642 is selectively phosphorylated in the presence of SRKR. 



INTRODUCTION 

The BK large-conductance calcium- and voltage- 
activated potassium channel (Kc a l.l, MaxiK) is gated by 
both membrane depolarization and intracellular Ca 2+ 
(Butler et al., 1993) . BK channels are tetramers of four 
pore-forming a subunits (Shen et al., 1994), encoded 
by the Kcnmal (Slol) gene (Butler et al., 1993; Meredith 
et al., 2004). BK a subunits are comprised of seven 
transmembrane domains (S0-S6) that form the chan- 
nel pore and are responsible for voltage sensing and a 
cytoplasmic C-terminal domain that binds Ca 2+ (Yuan 
etal., 2012). BK channels are expressed in a wide variety 
of excitable and nonexci table cells (Tseng-Crank et al., 
1994) where they have heterogeneous physiological 
functions. In neurons, BK channels contribute to the 
fast afterhyperpolarization of the action potential (AP) 
(Storm, 1987), AP duration and frequency (Storm, 
1987; Montgomery et al., 2013), and neurotransmitter 
release (Robitaille and Charlton, 1992; Hu et al, 2001). 
In smooth muscle, BK currents limit contraction and 
promote relaxation (Brayden and Nelson, 1992; Brenner 
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et al., 2000b). In nonexcitable cells, BK channels regu- 
late K + homeostasis (Rieg et al., 2007) and hormone se- 
cretion (Petersen and Maruyama, 1984) . Based on this 
wide array of cellular mechanisms, the BK channel has 
been implicated in the physiology of hearing (Fuchs 
and Evans, 1990; Pyott et al., 2007), circadian rhyth- 
micity (Meredith et al., 2006), cardiovascular function 
(Sausbier et al., 2005) , urination (Meredith et al., 2004) , 
and locomotor function (Sausbier et al., 2004; Imlach 
etal., 2010). 

Unlike many other K + channels, the BK channel tetra- 
mer is encoded by a single gene. Given this, physiologi- 
cal diversity of BK currents is achieved through several 
mechanisms including formation of macromolecular 
complexes with auxiliary subunits (Orio et al., 2002; 
Yan and Aldrich, 2012), posttranslational modification 
of a subunits (Tian et al., 2008; Yan et al., 2008), and 
alternative splicing (Shipston, 2001). Alternative splic- 
ing of Kcnmal transcripts has been proposed as a major 
mechanism driving the diversity of BK currents. Within 
the cytoplasmic C-terminal domain, four splice sites 
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have been identified. Site 1 is near the end of the first 
regulator of K + conductance (RCK) domain, containing 
either no insert or a four-amino acid alternate exon 
(SRKR) (Tseng-Crank et al., 1994) . Site 2 is the stress 
axis-regulated exon (STREX) site, found within the 
linker between RCK1 and RCK2. This site can contain 
either no insert, a three-amino acid exon (IYF), a 61— 
amino acid exon (STREX1), the IYF and 61-amino acid 
exons (STREX2), or a 28-amino acid exon (Saito et al., 
1997; Xie and McCobb, 1998; Chen et al., 2005). Site 3 is 
located near the high affinity Ca 2+ -binding site (the "Ca 2+ 
bowl") within RCK2 and can contain either no insert or a 
27-amino acid alternate exon (Ha et al., 2000). Lasdy, 
alternate exons at site 4 produce three distinct C termini 
("VEDEC," "VYR," or "QEERL") (Saito et al., 1997). 

Functionally, inclusion of alternate exons at each 
site produces effects on current properties that range in 
magnitude and mechanism, with the STREX site (site 2) 
being the most widely characterized. Inclusion of 
STREX1 or STREX2 enhances BK currents by speeding 
activation, slowing deactivation, and shifting the voltage 
dependence of activation more negative compared with 
insertless variants (Xie and McCobb, 1998; Chen et al., 
2005). Depending on the Ca 2+ concentration and the 
background splice variant sequence, inclusion of 
STREX can left-shift the V 1/2 up to 50 mV (Chen et al., 
2005) . Conversely, inclusion of the SRKR insert at site 1 
was found to reduce BK currents by either right-shifting 
the V1/2 (Tseng-Crank et al., 1994) or by reducing volt- 
age sensitivity (Rosenblatt et al., 1997). Inclusion of the 
Ca 2+ -bowl exon at site 3 does not have an appreciable 
effect on steady-state current properties; however, it 
does speed activation kinetics in a Ca 2+ -dependent man- 
ner (Ha et al., 2000) . Finally, the C-terminal exons so far 
appear to regulate trafficking and cell surface localiza- 
tion for the channel rather than significandy affecting 
current properties (Kim et al., 2008; Singh et al., 2013) . 

In addition to alternate splicing, posttranslational 
modifications, such as phosphorylation, have also been 
shown to underlie some of the physiological diversity in 
BK currents. The BK channel has ^200 potentially 
phosphorylated serine and threonine residues. How- 
ever, recent mass spectrometric analysis of rat brain BK 
protein detected only ^30 endogenously phosphory- 
lated residues (Yan et al., 2008). Some of these phos- 
phorylated residues were found within alternative exons 
at sites 2, 3, and 4 (Yan et al., 2008), suggesting that the 
combinatorial effects of exon-specific phosphorylation 
could further contribute to BK current properties. Con- 
sistent with this, inclusion of the STREX insert (site 2) 
changes the phosphorylation dependence of activation. 
Phosphorylation of S869 in channels lacking the STREX 
insert increases channel activity, whereas phosphoryla- 
tion of a serine within the STREX exon (S4 S trex) inhib- 
its channel activity (Tian et al., 2001). Although the 
phosphorylation dependence of STREX has been the 



most comprehensively studied, the functional impact 
of phosphorylation at other alternate exons has not 
been addressed. 

In this study, we cloned and characterized two novel 
Kcnmal variants (BKo and BK S rj<r) from the mouse su- 
prachiasmatic nucleus (SCN) of the hypothalamus, the 
brain's circadian clock. We chose to identify BK chan- 
nel sequences from SCN because (a) BK currents un- 
dergo dynamic transcriptional regulation in this tissue, 
suggesting that multiple splice variants may be expressed 
(Panda et al., 2002; Pitts et al., 2006); (b) distinct AP 
waveforms have been characterized in SCN neurons 
that are hypothesized to exert differential influence on 
intrinsic currents (Colwell, 2011; Montgomery and 
Meredith, 2012) ; and (c) SCN expresses BK modulatory 
(3 subunits (Montgomery and Meredith, 2012). The BKo 
and BKsrkr variants were identical except for the pres- 
ence of a four-amino acid insert (SRKR) at splice site 1 . 
To determine the functional effects of the SRKR insert, 
we expressed each variant with and without the (34 ac- 
cessory subunit in HEK293 cells, eliciting currents in 
response to standard voltage jumps and representative 
"day" and "night" SCN AP waveform commands. Next, 
using a combination of phosphatase treatment and site- 
directed mutagenesis, we addressed the role of two ser- 
ine residues, predicted to be phosphorylated in BK SR kr, 
in generating the distinct current properties produced 
by the BKo and BK S rkr variants. Lasdy, using SCN AP 
commands and expression of (34 to reconstitute activa- 
tion of BK currents under physiologically relevant con- 
ditions, we found that the differences in BK currents 
were expressed in response to the type of AP command, 
the splice variant expressed, the presence of the (34 sub- 
unit, and the phosphorylation status of S642. 

MATERIALS AND METHODS 

Cloning BK 0 and BK S r<r variants from SCN 

Total RNA was isolated from 10 pooled mouse SCNs harvested 
at zeitgeber time (ZT)6 (6 h after lights on), as described in 
Meredith etal. (2006). Total RNA was column purified (RNAeasy; 
QIAGEN), and cDNA was synthesized from 10 ug of total RNA 
(Superscript II; Invitrogen) using a gene-specific reverse primer 
(5'-GGTGACCATCATTCTCCTCAAAG-3'). PGR reactions ampli- 
fying the N-terminal ( 5 '-ATGGCTGTTGATGGGTGTTCGGG-3 ' 
and 5'-AGGCCCCCGAAGAAAGTCACCA-3') or a single product 
containing splice sites 1-4 ( 5 '-GAGTACAAGTCTGCCAACAG-3 ' 
and 5'-CATTCAAATCAAGCCCATGAGTACCC-3') were per- 
formed on 10 pi cDNA (Expand Taq; Roche). Products were gel 
purified, subcloned (pGEMTeasy; Promega), and sequenced to 
determine the presence of alternate exons at splice sites 1-4. Two 
PGR products were obtained with unique cDNA sequences, BKo 
and BKsrkr. 

A series of BK channel expression vectors containing the step- 
wise sequence additions and deletions of exons were introduced 
into mbr5 (Butler et al., 1993) in pcDNA3.1. In this study, mbr5 
contained an N-terminal myc tag and EYFP (716 bp) inserted in 
the RGK2 domain at nucleotide position (nt) 2032 (GenBank 
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accession no. KF530038; provided by R. Brenner, The University 
of Texas at San Antonio, San Antonio, TX). Construct 2 (C2; 
GenBank accession no. KF530044; Fig. SI) contained an inser- 
tion of an N-terminal myc tag followed by 195 bp of additional 
Kcnmal N-terminal sequence encoding the "MANG" alternate 
start site. The inserted sequence was synthesized as a minigene 
(Integrated DNA Technologies), digested with KpnI/Agel, and 
subcloned into mbr5/pcDNA3.1 at Kpnl and BssHII sites. 

For all subsequent constructs, splice site insertion or deletions 
were performed using site-directed PGR mutagenesis (Qui- 
kChange; Agilent Technologies) according to the manufacturer's 
protocol. For construct 3, three amino acids (IYF) were deleted 
from the STREX site by amplifying construct 2 with a primer con- 
taining a deletion of nts 2094-2102 (GenBank accession no. 
KF530036). For construct 4, 180 bp encoding the "VEDEC" C ter- 
minus was deleted by PGR from construct 3 (nts 4236-4415), and 
21 bp encoding the "RKEMVYR" C terminus was introduced in a 
second PGR amplification (GenBank accession no. KF530042). 
BKq was generated by the addition of 78 bp encoding the 27- 
amino acid Ca 2+ -bowl exon at nt 3564 in construct 4 (GenBank 
accession no. KF530040). BK S rkr was generated by the addition of 
12 bp encoding the SRKR exon at nt 1929 (GenBank accession 
no. KF530041). 

Quantitative PCR 

Total RNA was isolated from SCN tissue at 6 h after lights on 
(ZT6; "day") and 7 h after lights off (ZT19; "night") with the 
RNeasy mini kit (QIAGEN) as described previously (Montgomery 
and Meredith, 2012). Primers were as follows: SRKR, (F) 5'-GTTT- 
TGTGAAGCTTAAGCTCCTGAT-3 ' and (R) 5 '-AGAGCCGAAG- 
CCGAAAGC-3'; exon 16, (F) 5 '-GTTTTGTGAAGCTTAAGCTC- 
CTGAT-3' and 5 ' -TCTCTCTGTTGGC AGACTTGTAC-3 ' . Real- 
time qPGR data were analyzed using the 2~ AACq method (Livak 
and Schmittgen, 2001). SRKR and exon 16 Gq values for day were 
normalized to the night values. 

Predicted phosphorylation sites and mutagenesis 
of BKo and BKsrkr 

Sequence analysis for potential phosphorylation sites was per- 
formed online using NetPhos 2.0 (Blom et al., 1999). Predicted 
consensus phosphorylation sites received a score >0.9. 

Mutation of predicted phosphorylation sites S642 and S644 
was performed using site-directed PGR mutagenesis (QuikChange; 
Agilent Technologies) . BKo-S642D contained the mutation AGC to 
GAC at nts 1923-1925 (GenBank accession no. KF530040) . BKsrkr- 
S642A contained the mutation AGC to GCC at nts 1923-1935 
(GenBank accession no. KF530041). BK SRKR -S644A contained the 
mutation AGC to GCC at nts 1929-1931 (GenBank accession no. 
KF530041). All mutations were verified by sequencing. 

HEK cell patch-clamp electrophysiology 

HEK293T cells were transfected with BK expression constructs, 
with or without B4 in pcDNA3.1 (Wang et al., 2009). 1 pg 
DNA/35-mm dish was transfected for a subunits only, and 1 ug a 
plus 1.7 pg B4 was transfected for cotransfections, using Lipo- 
fectamine 2000 (Invitrogen) . Cells were plated on glass coverslips 
3-5 h later and recorded from 20 to 30 h after transfection. Patch- 
clamp recording was performed using the voltage-clamp mode in 
the inside-out patch configuration (Axopatch; Molecular De- 
vices) using thin-walled borosilicate pipettes with resistances of 
2-5 Mil. Recordings were made at room temperature, and the 
data were acquired at 50 kHz and filtered at 10 kHz. Total [Ca 2+ ]i 
was calculated using Webmax C Standard software and the appro- 
priate addition of CaCl 2 to achieve the indicated free Ca 2+ . For 
0 Ca 2+ solutions, 5 mM EGTA was used, where free Ca 2+ is ex- 
pected to be <1 nM under these conditions (Bao et al., 2002). For 
100-pM Ca 2+ solutions, no EGTA was used. For inside-out patch 



recordings, the pipette (extracellular) solution was composed of 
(mM): 140 KMeSO s , 2 KC1, 2 MgCl 2 , and 20 HEPES, pH 7.3. The 
bath (intracellular) solution was composed of (mM): 140 
KMeS0 3 , 2 KC1, 20 HEPES, and 5 HEDTAfor 1 and 10 pM of free 
Ca 2+ . To elicit BK currents using voltage-jump protocols, patches 
were stepped from holding potentials of — 100 to — 150 to +350 mV 
(in 10-mV increments) for 20 ms (oi alone) or 50 ms (a + B4), 
stepping back to —80 mV for 10 ms (a alone) or 20 ms (a + B4) to 
generate tail currents. Conductance (G) was calculated from tail 
currents, normalized to the highest conductance calculated for 
each patch, and plotted against membrane potential to generate 
G-V curves. Boltzmann functions were fitted to curves in Origin 8.5 
(OriginLab) using the equation: y = (Ai - A 2 )/(l + e (x ~ " 0 >/d " + A 2 . 
Single-exponential functions were fitted to current traces in 
pCLAMP (Molecular Devices) to determine the time constants of 
activation, deactivation, and inactivation, using the equation /(x) = 
A i e(-t/T i ) +C. 

For dephosphorylation experiments, calf intestinal alkaline 
phosphatase (New England Biolabs, Inc.) was diluted to 10 U/ml 
in the intracellular bath solution. After recording baseline cur- 
rents, phosphatase-containing solution (Alk P) was perfused into 
the bath, and posttreatment currents were recorded 5 min later. 

Activation of BK currents by AP commands 

AP command waveforms were originally recorded in whole-cell 
current-clamp mode from acute SCN slices during the day and 
night of the circadian cycle as described previously (Montgomery 
and Meredith, 2012). Day and night waveforms representative of 
the mean parameter values (see Results and Fig. 2 A) were used to 
elicit BK currents in whole-cell recording mode. The pipette (in- 
tracellular) solution was identical to the intracellular solution 
used in inside-out recordings, except with free Ca 2+ adjusted to 
50 pM. The bath (extracellular) solution consisted of (mM): 
145 NaCl, 5 KC1, 2 CaCl 2 , 1 MgCl 2 , and 10 HEPES, with pH ad- 
justed to 7.2 with NaOH. Cells were held at -100 mV. The ap- 
plied voltage stimulation protocol consisted of a depolarizing 
prepulse to +160 mV for 20 ms, followed by either day or night 
physiological AP commands that included holding the cells at the 
appropriate physiological resting membrane potentials (—48.6 mV 
for day AP commands and —45.0 mVfor night AP commands) for 
50 ms before the actual AP commands (Fig. S3 A). Series resis- 
tance was compensated up to 60%. 

For comparisons between splice variants, peak AP-elicited cur- 
rents were normalized to the +160-mV prepulse (Fig. S3 B). For 
AP-induced BK currents, peak current was the difference between 
the baseline current at —100 mV and the maximum current ob- 
served in response to the AP command. Current half-width was 
the duration at 50% of the maximal current magnitude, mea- 
sured between the current maximum and the afterhyperpolariza- 
tion current. Charge transfer was calculated as the area of the 
normalized current from baseline (at resting membrane poten- 
tial) to the maximum of the AP-induced current. 

Statistics 

Two-way ANOVAs with Bonferroni post-hoc analysis were used to 
determine statistical significance for differences between Vi/ 2 
values and day versus night AP-induced current parameters for 
the BK variants, as indicated in the figure legends. One-way 
ANOVA was used for comparisons between wild-type and mu- 
tant AP-elicited currents (Fig. 4 H) , and a t test was used to com- 
pare currents elicited by square voltage pulses in Fig. 4 G. 
Factorial ANOVAs with Bonferroni post-hoc analyses were used 
to determine statistical significance for T act and T deact between 
constructs elicited at different Ca 2+ across a range of voltages. 
Statistical significance was achieved if P < 0.05. Mean values in 
figures are given as ±SEM. 
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Online supplemental material 

Fig. SI demonstrates that the addition of the alternate N termi- 
nus, deletion of residues IYF at site 2, and the addition of the 
alternate C terminus (RKEMVYR) have little effect on BK cur- 
rent properties. Fig. S2 demonstrates the functional differences 
produced by the addition of the Ca 2+ -bowl exon at site 3 (BKo) 
and the addition of SRKR at site 1 (BK SRKR ) . Fig. S3 shows the 
prepulse voltage step used to normalize the AP-induced BK cur- 
rents in relation to the AP voltage command and the BK cur- 
rents elicited by each command to scale. Fig. S4 compares the 
amino acid sequences and secondary structures in the SRKR 
region from the two isoforms of the BK channel that have been 
crystallized to date. Table SI contains the Vi /2 values of all the 
constructs used in this study. The online supplemental mate- 
rial is available at http://www.jgp.org/cgi/content/full/jgp 
.201311072/DC1. 

RESULTS 

Identification of two novel BK channel splice variants 

To investigate BK current properties produced by 
native BK channels, we amplified full-length Kcnmal 
cDNAs from SCN, a tissue that exhibits dynamic expres- 
sion of Kcnmal transcripts (Panda et al., 2002; Pitts 
et al., 2006) . Two novel variants were identified (BKo and 
BK S rkr; Fig. 1 A) , containing an alternate exon combi- 
nation that differed at four splice sites from previously 
identified BK variants. Both variants contained an alter- 
nate translation start encoded by additional sequence 
in the extracellular N terminus of the channel ("MANG"; 
Fig. 1 A) . In addition, both variants lacked any insert at 
the STREX site (site 2), contained the 27-amino acid 
Ca 2+ -bowl insert (site 3), and contained the same alter- 
native C terminus, RKEMVYR (site 4) . However, the 
variants differed from each other at splice site 1 , located 
within RCK1 , where BK S rkr contains a four-amino acid 
insert, SRKR, that is absent in BKo. Stepwise addition of 
each alternate sequence into a BK channel lacking in- 
serts at splice sites 1 and 3, and replacements of exons 
at sites 2 and 4, revealed that the addition of the SRKR 
insert had the greatest impact on BK current properties 
(Figs. SI and S2). 

To determine whether the expression of BK transcripts 
containing the SRKR exon differed from the previously 
characterized nighttime-phased expression pattern of 
Kcnmal in SCN (Meredith et al., 2006; Montgomery 
and Meredith, 2012), exon 16 (the constitutive exon im- 
mediately preceding SRKR) and SRKR levels were deter- 
mined by quantitative RT-PCR from day and night 
SCNs (Fig. 1 B). The abundance of transcripts contain- 
ing exon 16 was greater at night compared with day 
(Fig. 1 B), similar to other constitutive Kcnmal exons 
(Panda et al., 2002; Pitts et al., 2006; not depicted). In 
contrast, the abundance of transcripts containing SRKR 
was higher during the day (Fig. 1 B), suggesting that 
SRKR-containing transcripts form a distinct type of BK 
channel that is regulated differentially with respect to 
circadian phase. 



SRKR alters kinetic and steady-state BK current properties 

Given the difference in abundance of SRKR-containing 
transcripts in day and night SCN tissue from the normal 
Kcnmal expression pattern, we looked at the functional 
impact of this exon on BK currents. BKo and BK SRKR were 
expressed in HEK293 cells, and macroscopic BK currents 
were recorded from inside-out patches in symmetrical 
K + at 0, 1, 10, and 100 pM of intracellular Ca 2+ . Expres- 
sion of both variants produced BK currents that acti- 
vated with increasing depolarization and Ca 2+ (Fig. 2 A) . 
Differences in current properties were apparent at 0, 1, 
and 100 pM, where the BK SRKR G-V relationship was 
right-shifted compared with BKo, but not at 10 pM Ca 2+ 
(Fig. 2 B). At 100 pM Ca 2+ , a condition where all of the 
high affinity Ca 2+ -binding sites are expected to be occu- 
pied, BK SRKR currents were harder to activate relative to 
BKo, with the Vi /2 right-shifted by 19 mV compared with 
BKq (Fig- 2 C) . This effect was also seen at 1 pM Ca 2+ 
(+35-mV right-shift), as well as in the absence of Ca 2+ 
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Figure 1. Native BK channel splice variants expressed in the 
SCN. (A) Schematic of the BK channel. Both BKo and BKsrkr 
contain an extended N-terminal region ("MANG"), no insert at 
site 2 (STREX), a 27-amino acid insert at site 3 ("Ca 2+ bowl"), 
and the "VYR" C terminus. In addition, BKsrkr includes the SRKR 
insert at splice site 1. (B) The day (open bars) and night (closed 
bars) expression of the alternate exon SRKR differs from the con- 
stitutive exon 16. Relative expression levels were determined by 
quantitative RT-PCR and normalized to the nighttime value for 
each reaction. Expression of transcripts containing exon 16 was 
significandy higher at night (n = 3; P = 0.0023; t test), whereas 
expression of SRKR-containing transcripts was higher during the 
day(n=3;P = 0.015;itest). 



588 



Functional effects of SRKR exon on BK currents 



(+49-mV right-shift) . There was no change in the slope 
factor (voltage sensitivity) between BKo and BK S rkr G-Vs 
at each Ca 2+ concentration. These data show that inclu- 
sion of the SRKR exon shifts the voltage dependence 
of activation to more positive values across a range of 
membrane potentials. Surprisingly, there was no signifi- 
cant difference in the V1/2 values of BK S rkr and BKo 
G-Vs at the intermediate Ca 2+ concentration of 10 uM 
(Fig. 2, Band C). 

There were also differences in the kinetics of activa- 
tion and deactivation between BKo and BK S rkr currents. 
At 0, 1, and 100 uM Ca 2+ , the time constants of activa- 
tion (t^) for BKsrkr currents were increased compared 
with BKo (Fig- 2 D). This slowing of activation was par- 
ticularly pronounced at 100 pM Ca 2+ at more negative 
voltages (Fig. 2 E). As with the steady-state G-V data, 
there was no difference in BKo and BK SRKR current acti- 
vation time constants in the intermediate Ca 2+ concen- 
tration of 10 pM. BKsrkr currents showed decreased 
Tdeact at 100, 1, and 0 pM Ca 2+ with no significant effect 



at 10-pM Ca 2+ concentrations (Fig. 2 F). Overall, at satu- 
rating (100 pM) Ca 2+ , the right-shifted V1/2, increased 
activation time constants, and decreased deactivation 
time constants would be expected to reduce native 
BKsrkr currents compared with BKo. The same effect is 
expected at low Ca 2+ (0 and 1 pM) , where BK S rkr cur- 
rents were similarly right-shifted and slower to activate 
compared with BKo, but absent at intermediate Ca 2+ 
(10 pM), where no significant differences in current 
properties were observed. 

SRKR is required for phosphoregulation of serine 642 

BKo and BK SR kr splice variants differ by only four amino 
acids, yet they have substantial differences in both 
steady-state and kinetic properties. Given the presence 
of a serine within the SRKR insert, we hypothesized that 
phosphorylation of this residue (S644; Fig. S4) could 
contribute to the difference in BKo and BK S rkr current 
properties. Analysis of SRKR and the surrounding 
sequence in BK S rkr with a phosphorylation prediction 



1 \iM Ca 2 - 



10 uM Ca 2+ 



Figure 2. Voltage-clamp re- 
cordings and analysis of BKo 
and BK SR kr currents from vari- 
ants expressed in HEK293 cells. 
(A) Current traces from BKo 
and BKsrkr patches elicited from 
voltage jumps at the indicated 
intracellular Ca 2+ . (B) Normal- 
ized G-V relationship for BKo 
and BKsrkr at different Ca 2+ . 
There was no significant differ- 
ences in the slopes of the G-V 
curves between BKo and BK SRK r 
at each Ca 2+ concentration (P = 
0.07; two-way AN OVA) . (C) Plot 
of V1/2 versus Ca 2+ exemplifying 
the difference in voltage de- 
pendence of activation between 
BKo an d BK SRK r. V 1/2 values were 
significantly different between 
BKo and BK SRKR at 0 Ca 2+ (198 vs. 
247 mV; P < 0.01; n = 8; two-way 
ANOVA), 1 pM Ca 2+ (169 vs. 
204 mV; P < 0.01; n = 8; two-way 
ANOVA) , and 1 00 pM Ca 2+ (-11 
vs. 8 mV; P < 0.01; n = 8; two-way 
ANOVA). No significant diffe- 
rence in Vj/2 was detected at 
10 pM Ca 2+ (64 vs. 67 mV; P = 0.1; 
n= 8; two-way ANOVA) . (DandE) 
Plot of T act versus voltage for BKo 
and BKsrkr at 0, 1, and 10 pM 
Ca 2+ (D) and 100 pM Ca 2+ (E). 
Activation of BKo was signifi- 
cantly faster than BKsrkr at 0, 1, 
and 100 pM (P < 0.01; factorial 
ANOVA), but not at 10 pM Ca 2+ , 
where no differences in kinetics 

were observed. All symbols as in B. (F) Plot of T deact versus voltage for BKo and BKsrkr at 0, 1, 10, and 100 pM Ca 2+ . Deactivation of BKo 
was significantly slower than BK SR kr with 0, 1, and 100 pM Ca 2+ , but not with 10 pM Ca 2+ (P < 0.01; factorial ANOVA). All symbols as in 
B. In cases where error bars are not present, the bars were smaller than the data symbol. 
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algorithm (Blom et al., 1999) indicated a high probability 
for S644 phosphorylation (NetPhos score of 0.991). S642, 
a constitutive residue preceding the SRKR insert (Fig. S4) , 
also had a high probability score for phosphorylation 
(NetPhos score of 0.994) . In contrast, in the BKo sequence, 
which lacked SRKR, the constitutive S642 had a signifi- 
candy lower probability of phosphorylation (0.557) . These 
results predicted that phosphorylation of S642 and/ or 
S644 may occur in BK S rkr, but not in BKo, channels. 

First, to determine if net phosphorylation of BK S rkr 
could explain the differences in current properties com- 
pared with BKq, we applied alkaline phosphatase to the 



intracellular side of macropatches. At 100 pM Ca , 
the BKsrkr G-V curve was right-shifted compared with 
BKo. In the presence of phosphatase, the BK S rkr Vi/g 
was left-shifted by 14 mV, making it similar to that of 
BKo and significantly different from that of BK S rkr (Fig. 3, 
A and B) . This result suggested that dephosphorylation 
increased BK S rkr channel activity to match that of BKo. 
In contrast, the application of phosphatase had no ef- 
fect on BKo currents (Fig. 3, A and B), confirming that 
the SRKR insert is required to render the currents sus- 
ceptible to phosphatase. These experimental results are 
consistent with the high phosphorylation probability 




K)0~ 



Too" 



E 14 
12- 
10- 

1 * 
J* 

4 

2\ 



'SRKR 

BK +Alk P 



200 



150 



S100 



50 



-50 



50 100 



-50 



Membrane Potential (mV) 




200- 



150- 



> 

BK SRKR -S644A 
BK SRKR -S642A 
BK„-S642D 



S100- 



-50 



50 



S644A 

Too 



50- 



-50 



Membrane Potential (mV) 



BK^-S642A 




F 3.5 

3.0 
_2.5 
4.2.0 
h'1.5 

1.0 



0.5' 



3 °M 



0 



50 



embrane Potential (mV) 



100 




10 



Ca 2 



Too 




10 



Ca 2 



'100 




^2~5o ' ~aEq ' ^Too 

Membrane Potential (mV) 



-50 



Figure 3. BKo and BKsrkr current 
properties in patches treated with phos- 
phatase or from phospho-site mutants. 
(A) G-V relationship of BKo and BKsrkr 
at 100 uM Ga 2+ before and after the 
application of calf intestinal alkaline 
phosphatase (Alk P). Alk P caused a 
leftward shift of the BK SRKR G-V but not 
of BKo G-V. (B) Plot of Vi/ 2 versus Ca 2+ , 
illustrating that dephosphorylation of 
BKsrkr produces a BKo-like phenotype 
in the voltage dependence of activa- 
tion. After treatment with phosphatase, 
the Vi /2 was significantly left-shifted 
with BKsrkr at 1 uM Ca 2+ (204-155 mV; 
p = 0.01; n = 8; two-way ANOVA) and 
at 100 uM Ca 2+ (8 to -6 mV; P = 0.03; 
n = 8; two-way ANOVA), but not with 
BKo at either 1 pM Ca 2+ (169-156 mV; 
P = 0.25; n = 8; two-way ANOVA) or 
100 pM Ca 2+ (-11 vs. -7 mV; P = 0.53; 
n = 8; two-way ANOVA), making the 
phosphatase-treated BK SRK r not sig- 
nificandy different from BKo (P = 0.7; 
n = 8). Symbols are as in A. (C) G-V 
relationship at 1 00 pM Ca 2+ of specific 
phospho-mutants used to determine 
the site of differential phospho-regula- 
tion of BKsrkr. Phospho-mutant BK SRK r 
S644A, the S in SRKR, had no effect 
on G-V relationship, whereas phospho- 
mimetic BK SRK r S642D, found within 
a constitutive portion of the channel, 
was found to be necessary and suffi- 
cient to induce that change in current 
properties. (D) Plot of V1/2 versus Ca 2+ 
for BK constructs shown in C. At 1 and 
100 pM Ca 2+ , the V V2 values cluster 
into two groups, with those of BKsrkr, 
BKsrkr-S644A, and BKo-S642D further 
At 1 pM Ca 2+ , there were no significant 
S— 12; two-way ANOVA) , between BKsrkr- 



right-shifted than those of BKo, BK SR kr-S642A, and BK SR kr-S642A/S644A. Symbols are as in C 
differences in the V V2 values between BKsrkr and BKsrkr-S644A (204 and 206 mV; P = 0.29; n = 
S642A and BKo (166 and 169 mV; P = 0.29; n = 8-10; two-way ANOVA) , between BKo-S642D and BKsrkr (209 vs. 204 mV; P = 0.76; n = 8; 
two-way ANOVA), or between BK SR kr-S642A/S644A and BKo (167 vs. 169 mV; P = 0.25; n = 8-15; two-way ANOVA). At 100 pM Ca 2+ , 
there were no significant differences in the V1/2 values between BKsrkr and BKsrkr-S644A (6 and 8 mV; P = 0.3; n = 8; two-way ANOVA) , 
between BK SRK r-S642A and BKo (V V2 values both - 1 1 mV; P = 0.6; n = 8; two-way ANOVA) , between BKo-S642D and BKsrkr ( 1 1 vs. 8 mV; 
p = 0.4; n = 8; two-way ANOVA) , or between BK SRK r-S642A/ S644A and BKo ( - 15 vs. - 1 1 mV; P = 0.6; n = 8; two-way ANOVA) . (E and F) 
Plots of T act (E) and T deact (F) versus voltage for BKo, BKsrkr, and BKsrkr-S642A at 100 pM Ga 2+ . Activation kinetics of BK SRK r-S642A were 
not significantly different to those of BKo (P = 0.34; n = 8) but were significantly slowed compared with BKsrkr (P = 0.02; n = 8-9; facto- 
rial ANOVA) . Similarly, deactivation kinetics of BK SRK r-S642A were not significantly different to those of BKo (P = 0.49; n = 8; factorial 
ANOVA) but were significantly slowed compared with BK SR kr (P < 0.01; n = 8-9; factorial ANOVA) . 
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score for the BKsrkr variant and low probability score 
for the BKo variant. 

Next, to determine if a serine residue was required 
within the SRKR insert, we mutated S644 to an alanine. 
We hypothesized that currents produced from BK SRK r- 
S644A mutants would resemble both BKo and the 
phosphatase-treated BK SRK r currents. Surprisingly, 
BK S rkr-S644A showed no difference in the G-V relation- 
ship compared with wild-type BK S rkr (Fig- 3, C and D) . 
Furthermore, phosphatase treatment of BK SR kr-S644A 
still left-shifted the G-V curve to match that of BKo (not 
depicted) . Collectively, these results suggested that de- 
phosphorylation is occurring at a residue other than 
S644 to produce the observed functional effects. 

Next, we focused on S642, a constitutive serine pres- 
ent in both BKo and BK SRK r variants. We hypothesized 
that in the absence of SRKR, S642 would be resistant to 
phosphorylation. Conversely, the presence of the SRKR 
insert may render S642 permissive to phosphorylation, 
leading to the observed right-shifted G-V curve for 
BKsrkr currents. To test this prediction, S642 was mu- 
tated to alanine. Confirming the role of this serine in 
the functional current differences, the G-V curve of 
BKsrkr-S642A matched that of BKo (Fig. 3, C and D). 
Furthermore, the introduction of the converse phos- 
phomimetic mutation into BKo, BKo-S642D, produced a 
right-shifted G-V curve that matched that of BK SRK r 
(Fig. 3, C and D) . These results identify S642 as the criti- 
cal residue underlying the differences in BK S rkr current 
properties compared with BKo and suggest that phos- 
phorylation of S642 decreases channel activity. 

Mutating both serine residues to alanine in BK SRK r 
(BK SRK R-S642A/S644A) resulted in G-V curves that 
matched those of BKo (Fig. 3, C and D), and no further 
left-shift of the G-V was produced by phosphatase treat- 
ment of this mutant (not depicted) . These data confirm 
that no additional phosphorylation sites (sensitive to 
this phosphatase) were involved in the G-V differences 
between BKo and BK SRK r. Quantitatively similar results 
were obtained when the phosphatase and mutagenesis 
experiments were repeated at 1 uM Ca 2+ (Fig. 3 D). 
However, at 10 pM Ca 2+ , in which no functional differ- 
ences in the currents produced by BKo and BK S rkr were 
observed (Fig. 2), neither phosphatase treatment nor 
mutagenesis altered the G-V curves (Fig. 3, B and D). 

In addition to steady-state current properties, muta- 
tion of S642 also impacted current kinetics. We found 
that the time constants of both activation and deactiva- 
tion of BK S rkr-S642A were similar to those of BKo (Fig. 3, 
E and F) . Thus, both the steady-state and kinetics prop- 
erties of BKsrkr lacking the S642 phosphorylation site 
resemble those of BKo. 

SRKR reduces BK currents elicited by AP commands 

Much of the previously published work on BK current 
properties produced by specific variants has examined 



currents in response to standard square voltage-pulse 
protocols. Few studies have addressed the properties of 
BK currents elicited by a more physiologically relevant 
stimulus, such as an AP. Therefore, we recorded BK cur- 
rents elicited from two distinct types of AP waveforms 
(day and night) , based on recordings from SCN neu- 
rons where BKo and BK SRK r are endogenously expressed. 
These experiments examined whether splice variant- 
based differences in BK currents were still apparent in 
response to physiologically relevant activation. We also 
examined currents elicited under these conditions co- 
expressed with (34, the major neuronal auxiliary sub- 
unit, which is known to be expressed within the SCN 
(Montgomery and Meredith, 2012). 

The representative day and night AP waveforms used 
as a voltage command for BK current activation were 
obtained from whole-cell recordings of spontaneously 
firing neurons during the peak (day) and trough 
(night) of neural activity (Montgomery and Meredith, 
2012) . Day and night AP waveforms differed significandy 
in several parameters (Montgomery and Meredith, 
2012), including peak depolarization (14.9 mV for day 
compared with —2.2 mV for night), half-width (2.1 vs. 
3.0 ms, respectively), and baseline membrane potential 
(—48.6 vs. —45.0 mV, respectively). These distinct AP 
properties were expected to exert a differential drive on 
BK activation, conditions that might alter the expres- 
sion of the differences observed under square voltage- 
jump conditions between BKo and BK S rkr currents. 

Day and night AP command waveforms (Fig. 4 A) were 
used to elicit currents in physiological ion concentra- 
tions (low internal K + and high external Na + ). Fig. 4 B 
shows AP-elicited currents from BKo and BK S rkr, elicited 
from day or night AP command waveforms at 50 pM 
Ca 2+ , mimicking the high local Ca 2+ concentration expe- 
rienced by BK channels in vivo coupled with voltage- 
gated Ca 2+ channels (Fakler and Adelman, 2008) . In this 
configuration, the time course of the BK currents gener- 
ally followed the depolarizing phase of the AP command 
waveforms, mirroring the differences seen between the 
day and night AP properties. To directly compare BKo 
and BKsrkr properties across cells expressing different 
levels of each channel variant, the AP-elicited currents 
were normalized to a maximally activating depolarizing 
prepulse step (Fig. S3) . The mean AP-induced peak cur- 
rents from both variants were very low, ranging from 2 to 
8% of maximal current, depending on the variant and 
AP command waveform. Both BKo and BK S rkr currents 
elicited by day AP commands had significandy increased 
peak current amplitudes and decreased half-widths com- 
pared with currents elicited by night AP commands 
(Fig. 4, D and F) . There was no significant difference in 
charge transfer elicited by day and night AP commands 
(Fig. 4, E and F) . These data show that the day and night 
AP command waveforms produced differing drives on 
BK current activation under standardized conditions. 
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Next, we compared the properties of BKo and BK S rkr 
currents to each other, using either day or night AP 
commands (Fig. 4 C) . BKo and BK S rrr currents elicited 
by the same stimulus, either the day or night AP com- 
mand waveform, showed different properties (Fig. 4 C) . 
BK SRK r had significandy decreased peak currents (Fig. 4 D) , 
decreased charge transfer (Fig. 4 E), and increased cur- 
rent half-width (Fig. 4 F) compared with BKo. These re- 
sults are consistent with our findings using standard 
voltage jumps (Fig. 2 B) .The reduced BK S rkr currents 
compared with BKo elicited by AP commands could be 
caused by differences in the steady-state properties of 
the channels, kinetic properties, or a combination of 
both. To determine whether reduced steady-state acti- 
vation could explain the differences between BKo and 
BKsrkr AP-induced currents, we used a square pulse of 



40 ms to depolarize cells to +14.9 mV, the peak voltage 
achieved using the day AP command. We found no sig- 
nificant difference between BKo and BK SR kr currents 
under these conditions (Fig. 4 G), indicating that differ- 
ences in steady-state parameters are not responsible for 
the distinct AP-induced BKo and BK S rkr current proper- 
ties. This suggests that kinetic differences account for 
the major differences in currents produced by the two 
variants elicited by AP commands. 

Because S642 is critical for expression of the differ- 
ences between BKo and BK SRK r currents under standard 
voltage-jump protocols, we next sought to determine 
whether removal of the phosphorylation site S642 in 
the BKsrkr variant would render channels "BKrlike" in 
their responses to AP command waveforms. BK SRK r- 
S642A currents elicited by the day AP command were 
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Figure 4. BKsrkr and BKo cur- 
rents elicited from AP com- 
mands in HEK293 cells. (A) Day 
(red) and night (black) AP 
commands, recorded from SCN 
neurons. (B) Currents elicited 
from BKo and BKsrkr in physi- 
ological K + at 50 uM Ca 2+ in 
response to either day (red) or 
night (black) AP commands in A. 
(C) BKo (black) and BKsrkr (gray) 
currents elicited from day (left) 
and night (right) AP commands, 
normalized to the peak steady- 
state current in response to a 
+ 160-mV step. (D-F) Mean 
normalized peak current (D), 
charge transfer (E), and half- 
width (F) values for BKo and 
BKsrkr currents elicited by day 
(open bars) and night (closed 
bars) AP commands. For both 
splice variants, currents induced 
by day versus night AP com- 
mands had significantly greater 
peak current amplitude (P = 
0.04; two-way ANOVA) and sig- 
nificantly reduced half-width 
(P < 0.01; two-way ANOVA) 
compared with currents elicited 
from night APs. There was no 
significant difference in charge 
transfer between day and night. 
Comparing the values between 
splice variants, BK SRK r currents 
were reduced compared with 
BKo currents using both day and 
night AP commands, as shown 
by the increased peak current, 
charge transfer, and half-width 
(P < 0.01 for all BKo and BKsrkr 

comparisons; n = 9-11 for all conditions; two-way ANOVA). (G) Normalized peak currents for BKo and BK SRK r elicited from day AP 
commands were significantly different, whereas those elicited from a square voltage pulse of the same magnitude of the day AP were 
not significantly different (P = 0.4; n = 3-11; t test). (H) Normalized peak currents in response to day AP commands from BKsrkr-S642A 
were not significantly different from BKo (P > 0.5; one-way ANOVA) but were significantly different from BK SRK r (P = 0.05; n = 9-11; 
one-way ANOVA) . 
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significantly larger than wild-type BK S rkr currents, and 
were similar to BKo peak currents (Fig. 4 H) . This result 
suggests that phosphorylation of S642 in BK S rkr channels 
would be expected to alter BK currents activated during 
the AP, opening the possibility that the SRKR splice site 
could play a role in influencing neural activity. 

(34 enhances the functional differences between BK 0 
and BK SRK r 

In vivo, BK current properties are modified by associa- 
tion with accessory (3 subunits. In the SCN (Montgom- 
ery and Meredith, 2012) and the brain as a whole, (34 
is the predominant BK modulatory subunit (Behrens 
et al., 2000; Brenner et al., 2000a). Coexpression of (34 
with the mbr5 BK channel has been shown to exert 
complex effects, decreasing Ca 2+ sensitivity at low Ca 2+ 
and increasing Ca 2+ sensitivity at high Ca 2+ , with litde 
effect at 10 pM Ca 2+ (Behrens et al, 2000; Brenner 
et al., 2000a; Lippiat et al., 2003; Wang et al., 2006). (34 
also influences channel kinetics, slowing rates of both 



activation and deactivation (Behrens at al., 2000; 
Brenner et al., 2000a). To examine the influence of (34 
on each splice variant, and to determine whether (34 af- 
fects the functional differences observed between the 
splice variants, we coexpressed (34 with either BKo or 
BKsrkr. First, currents were recorded from inside-out 
patches in symmetrical K + solutions with standard 
square voltage jumps (Fig. 5 A), and G-V relationships 
and activation and deactivation kinetics were examined 
(Fig. 5, B-G). Similar to the previously reported ef- 
fects on other BK channel variants, coexpression of (34 
slowed activation of both BKo and BK S rkr at both 10 and 
100 pM Ca 2+ (Fig. 5, A and B) . At 10 pM Ca 2+ , (34 slowed 
the activation time constants of BKo and BK S rkr to a 
comparable degree. Thus, the activation time constants 
at 10 pM Ca 2+ were similar between the variants with 
and without (34 (Figs. 2 D and 5 B) . 

In contrast, at 100 pM Ca 2+ , although (34 slowed the 
activation of both BKo and BK S rkr, the relative differ- 
ences in activation kinetics between the variants in the 
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Figure 5. BKo an d BK SRKR currents from variants co- 
expressed with [54. (A) BKo and BK SR kr currents were elic- 
ited with voltage steps from -100 to 200 mV with 10 pM 
Ca 2+ . Current traces were normalized to I max elicited from 
either the a subunit alone (black) or coexpressed with (34 
(gray traces). (B and C) Plot of T act (B) or T deact (C) versus 
voltage for BKo an< f BKsrkr coexpressed with [54 at 10 and 
100 pM Ca 2+ . Symbols are as in C. Although no difference 
is found in kinetics between the two variants at 10 pM Ca 2+ , 
BK SRK R was significantly (P < 0.01 ; factorial ANOVA) slower 
to activate and faster to deactivate (P < 0.01) when coex- 
pressed with (34 at 100 pM Ca 2+ . (D and E) G-V relation- 
ships of BK, an <f BK SRKR coexpressed with (34 at 10 pM Ca 2+ 
(D) and 100 pM Ca 2+ (E). Vi/ 2 values were significantly dif- 
ferent between BKo and BKsrkr at both 10 pM Ca 2+ (151 vs. 
94 mV; P < 0.01; n=8; two-way ANOVA) and at 100 pM Ca 2+ 
(15 vs. -23 mV; P < 0.01; n = 8; two-way ANOVA). (F) Plot 
of AV1/2 versus Ca 2+ for each splice variant caused by (34 
expression (Vi/ 2 [a|3] — Vi /2 [c«]). Increasing Ca 2+ reduces 
the V1/2 shift caused by (34 for both variants, albeit to differ- 
ing extents. (G) Plot of AV1/2 versus Ca 2+ between variants 
(Vi/ 2 [BK SR kr] — V 1/2 [BK,]). The addition of (34 increases 
the difference in Vi/ 2 between splice variants. 



Shelley etal. 593 



presence of (34 were increased at positive voltages 
(Fig. 5 B) compared with activation kinetics in the ab- 
sence of (34 (Fig. 2 E). The net effect was a larger diver- 
gence in activation time constants between the variants 
when (34 was coexpressed with 100 pM Ca 2+ . 

(34 also slowed deactivadon of both splice variants at 
100 pM Ca 2+ , while sdll maintaining the faster deactiva- 
tion of BKsrkr compared with BKo at more depolarized 
potentials (Fig. 5 C) . In the absence of (34, differences 
in BKo and BK SRKR deactivation were most apparent at 
potentials hyperpolarized to approximately —100 mV 
(Fig. 2 B). However, in the presence of (34, increasing 
voltage caused a widening divergence of deactivation 
time constants between BKo and BK SRK r (Fig. 5 C). At 
10 pM Ca 2+ , (34 also slowed deactivation of both variants 
approximately fourfold. As with activation time con- 
stants, the deactivation time constants of the two vari- 
ants at 10 pM Ca 2+ in the absence of (34 were similar 
(Fig. 2 F), and (34 coexpression did not introduce any 
further relative differences between BKo and BK S rkr. 

Although (34 coexpression did not affect the differ- 
ences in kinetics of activation or deactivation between 
the variants at 10 pM Ca 2+ , it did influence the G-V rela- 
tionship. At 10 pM Ca 2+ , the presence of (34 right- 
shifted the Vi/9 values for BKo and BK S rkr by 30 and 
84 mV, respectively (Fig. 5, D and F, and Table SI) . Fur- 
thermore, although in the absence of (34 at 10 pM Ca 2+ 
no difference was seen between the G-V curves of BKo 
and BKsrkr (Fig- 2, B and E), coexpression with (34 now 
produced a 57-mV difference in V1/2 values (Fig. 5 G). 

In contrast to 10 pM Ca 2+ , at 100 pM Ca 2+ , (34 co- 
expression was found to have much reduced, but op- 
posing, effects on the G-V curves. For BKo, the addition 
of (34 left-shifted the V1/2 by 12 mV; however, in the case 
of BKsrkr, (34 caused a 7-mV right-shift (Fig. 5, E and F, 
and Table SI ) . Coexpression of (34 increased the differ- 
ences in V1/2 between the variants, although to a lesser 
extent than at 10 pM Ca 2+ (Fig. 5 G) . These results dem- 
onstrate that (34 coexpression acts to increase the func- 
tional differences between the steady-state properties of 
the splice variants at 10 and 100 pM Ca 2+ . 



We next examined the consequence of (34 coexpres- 
sion on the splice variant responses to day and night SCN 
AP waveforms. Whole-cell currents were recorded with 
physiological solutions (low internal K + and high external 
Na + ) at 50 pM of internal Ca 2+ . In response to night AP 
commands, (34 coexpression led to a significant twofold 
increase in BKo peak current compared with BKo alone 
(Fig. 6 B) . This effect was consistent with the left-shifted 
G-V from BKo/ (34 currents at a similar high Ca 2+ (Fig. 5 E) . 
In contrast, BK S rkr peak currents elicited by night AP 
commands were not significantiy enhanced by coexpres- 
sion of (34, agreeing with the modestly right-shifted G-V 
for BKsrkr/ (34. In addition, (34 had no significant effect on 
either BKo or BK S rkr currents elicited by day AP com- 
mands (Fig. 6, A and B) . These results show that (34 dif- 
ferentially affects the BKo and BK S rkr splice variants under 
distinct stimulus conditions, an effect that is consistent 
with results obtained in symmetrical K + conditions (Fig. 5) . 

Because BK SR kr peak currents were reduced compared 
with BKo in response to both day and night AP commands 
without (34, we compared the magnitude of the decrease 
in BKsrkr peak current in the presence of (34 (Fig. 6) . In 
response to the night AP command without (34, BK SRK r 
currents were decreased by 2.5-fold compared with BKo. 
In the presence of (34, this difference widened to three- 
fold (Fig. 6 B) . In contrast, the day AP command pro- 
duced a significant 2.8-fold reduction in BK S rkr current 
compared with BKo in the absence of (34, which was 
blunted to a 2.1-fold reduction by the addition of (34 
(Fig. 6 A). Furthermore, in the presence of (34, BK SRK r 
peak currents were no longer statistically different to BKo 
because of an increase in the variability with (34 coexpres- 
sion. Collectively, these results suggest that (34 plays a 
greater role in modulating the differences between BKo 
and BKsrkr currents in response to night AP commands 
than in response to day AP commands. 

DISCUSSION 

In this study, we identified two BK channel splice vari- 
ants containing novel combinations of alternate exons 
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Figure 6. (34 effects on AP-induced currents. (A and B) 
Normalized peak currents for BKo and BK SRK r elic- 
ited by day (A) and night (B) AP commands, ex- 
pressed with and without (34. In the absence of (34, 
BKsrkr currents were significandy reduced compared 
with BKq f° r both day and night AP commands (P < 
0.05; data replotted from Fig. 4 to facilitate compari- 
sons; two-way ANOVA). Coexpression of (34 had no 
significant effect on peak current from either splice 
variant elicited by day AP commands. In contrast, (34 
significandy increased BKo peak currents (P < 0.01; 
two-way ANOVA) , preserving the significant differ- 
ence in BKo + (34 and BKsrkr + (34 peak currents (P < 
0.01; n = 5—11; two-way ANOVA) . For clarity, signifi- 
cant differences are only marked for pairwise com- 
parisons involving (34. 
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(BKo and BK S rkr) , which differed from each other by a 
four-amino acid insert (Fig. 1) . We determined that the 
functional differences between BKo and BK SR kr currents 
were caused by the SRKR exon (Fig. 2) . Inclusion of the 
SRKR exon (BK S rkr) resulted in a large reduction of BK 
channel activity, demonstrated by a rightward shift in 
G-V relationships, slower activation rates, and faster de- 
activation rates. These effects were abrogated by dephos- 
phorylation only when SRKR was present and, conversely, 
could be reproduced in the absence of SRKR by a phos- 
phomimetic mutation (S642D) (Fig. 3). Overall, the 
changes in current properties caused by SRKR were 
larger than those caused by the inclusion of other al- 
ternate exons at sites 2, 3, and 4 (Figs. SI and S2, and 
Table SI) . When BK S rkr was examined in a reconstituted 
physiological system activated by neuronal AP com- 
mands in saturating Ca 2+ , we found that even with the 
slower rate of depolarization and the much shorter du- 
ration of AP commands compared with voltage jumps, 
BKsrkr still had reduced activity compared with BKo. This 
inhibition was manifested as a reduction in peak current, 
charge transfer, and current half-width during the AP 
stimulus compared with BKo and was eliminated by muta- 
tion of S642 (Fig. 4) . Collectively, the central finding of 
this study establishes the context-dependent phosphoryla- 
tion of S642 in the presence of the SRKR exon as a novel 
means for reducing BK activation during an AP stimulus. 

The functional significance of SRKR inclusion in BK 
transcripts parallels the known physiological role for BK 
currents in regulating neural activity in the SCN. BK cur- 
rents primarily regulate firing frequency at night in SCN 
neurons, correlating with the window of higher total 
BK channel expression and larger BK currents in wild- 
type SCNs (Panda et al., 2002; Meredith et al., 2006; 
Montgomery et al., 2013) . Yet BK currents have little im- 
pact on firing rate during the day, despite the presence 
of a baseline BK current during this time in wild-type 
SCN neurons (Meredith et al., 2006; Pitts et al., 2006). 
Our data suggest that part of the reduced influence for 
BK currents on SCN neuronal activity during the day could 
be mediated by an increase in SRKR-containing tran- 
scripts, leading to decreased AP-induced BK currents. 

To further understand how the changes in BK SRK r 
current properties lead to a reduction in AP-induced 
current, we considered the effect of SRKR on both ki- 
netic and steady-state parameters, finding alterations in 
both compared with BK,. Mutation of S642 to alanine 
(BK SRK r-S642A) converted BKsrkr steady-state and ki- 
netic properties to be more "BKo-like" (Fig. 3, C-F). 
However, the reduction in BK S rkr current compared 
with BKo was eliminated when a step command to the 
same peak voltage applied during the AP command was 
used (Fig. 4 G), suggesting that the reduced BK SRK r 
current elicited by APs may be largely based on kinetic 
differences between the variants. Consistent with this, 
BKsrkr currents exhibited slowed activation and faster 



deactivation at high Ca" + (Fig. 2), factors that could be 
expected to prominently reduce BK S rkr current activa- 
tion in response to a brief AP stimulus compared with 
BKo. Interestingly, when (34 was coexpressed, which fur- 
ther slows both activation and deactivation of BK SRK r, 
only the night AP recapitulated the significant differ- 
ence in peak current between the variants (Fig. 6) . In 
this case, the slower rise time of the night AP compared 
with the day AP (Fig. 4 A) may facilitate expression of 
the intrinsic difference between BKo and BK SR kr- Col- 
lectively, these data strongly support the idea that the 
expression of BKo and BK S rkr is expected to impact ex- 
citability in vivo. 

Variants with and without the SRKR insert have been 
previously cloned from human, chick, and turtie tissues 
(Tseng-Crank et al., 1994; Rosenblatt et al., 1997; Jones 
et al., 1999). The suppressive effect of SRKR on BK cur- 
rents in this study is consistent with some data from 
more limited investigations of the SRKR insert. In the 
human BK variants, SRKR right-shifted the V1/2 of 
-10 mV at 2.4, 75, and 405 uM Ca 2+ (Tseng-Crank et al., 
1994) . In another study, BK variants with and without 
SRKR isolated from the chick cochlea showed no sig- 
nificant differences in V1/2 values, but the SRKR-con- 
taining clone had reduced voltage sensitivity, indicated 
by a shallower G-V curve (Rosenblatt et al., 1997). In 
turtle cochlear variants, the presence of SRKR could 
either speed or slow activation and left- or right-shift 
G-V curves, depending on the exons at other splice sites 
(Jones et al., 1999). Comparing these data to our re- 
sults, we observed larger V1/2 right-shifts of 20-40 mV 
and no effect on the G-V slope. These differences may 
be attributable to the novel combination of exons con- 
tained within BKo and BK SR kr splice variants that may 
interact to modulate the functional effects of SRKR. For 
example, the human variant also contained the STREX 
insert, the chick variant contained an additional eight- 
residue insert in the RCK1 domain, and both lacked the 
27-residue Ca 2+ -bowl insert at splice site 3 (Tseng-Crank 
et al., 1994; Rosenblatt et al., 1997). In our study, the 
addition of the Ca 2+ -bowl exon also right-shifted the 
V1/2 (Fig. S2). This exon is present in both BKo and 
BKsrkr variants and may enhance the right-shifting ef- 
fect of SRKR in BK SR kr. Such functional interactions 
between alternative exons have been described previ- 
ously in BK channels cloned from Caenorhabditis elegans 
(Johnson et al., 2011). For example, insertions at a 
splice site within the RCK1 domain can modulate Ca 2+ 
sensitivity and activation kinetics, but only when inser- 
tions are present at further splice sites within the RCK1- 
RCK2 linker (Johnson etal., 2011). 

Mechanistically, how does inclusion of the SRKR 
exon regulate phosphorylation? Consistent with our re- 
sults showing that the dephosphorylation-dependent 
functional differences at S642 cannot be detected with- 
out SRKR present, insertion of SRKR likely alters the 
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structure of the local protein environment to render 
S642 susceptible to phosphorylation. Examination of 
the available crystal structures of BK channels supports 
this idea (Fig. S4). The structure of the cytosolic do- 
main of a human BK channel variant (Wu et al., 2010) 
contains the homologous residue to S642 but lacks the 
SRKR insert. In this structure, S642 is positioned at the 
extreme N terminus of a (3 strand immediately after a 
short five-residue loop toward the C-terminal end of the 
RCK1 domain. The only other isoform for which there 
is a crystal structure is a zebrafish isoform that contains 
the SRKR insert (Yuan et al., 2012). In this structure, 
the atomic coordinates corresponding to S642 and sev- 
eral surrounding residues, including the SR portion of 
SRKR, were not obtained. However, the structure sug- 
gests that S642 resides within an inter (3-strand loop; 
therefore, the presence of SRKR shifts S642 from within 
a (3-strand region to a loop (Fig. S4) . The shift of S642 
into this loop could better position it for accessibility to 
kinases for phosphorylation. In addition, SRKR inser- 
tion would distance S642 from the subsequent hydro- 
phobic residues ILI, the proximity of which reduces the 
predicted phosphorylation score in NetPhos 2.0. Inter- 
estingly, of the 30 phospho-residues detected in rat 
brain lysates (Yan et al., 2008), no peptide fragment 
containing S642 was detected at all (phosphorylated or 
not), presumably because of its small size, leaving open 
the investigation for BK SR kr variants phosphorylated at 
S642 in native tissues. 

How might SRKR and a phosphate group at S642 act 
to reduce channel activity? In other studies, phosphory- 
lation has been shown to either increase or decrease BK 
channel activity, depending on the kinase, residue posi- 
tion, and channel variant (Reinhart et al., 1991; Ling 
et al., 2000; Tian et al., 2001; Yan et al., 2008). In the 
case of S642 and SRKR, these residues are not located in 
the vicinity of any of the divalent ion-binding sites, nor 
are they close to the transmembrane regions that form 
the voltage sensors or the pore. In the absence of Ca 2+ , 
the rightward V1/2 shift of BKsrkr compared with BKo 
indicates that phosphorylation of S642 reduces BK 
channel activation via steps that are Ca 2+ independent, 
for example, those involved in intrinsic pore opening 
or the coupling of voltage-sensor movement to pore 
opening. However, direct effects on voltage-sensor 
movement are unlikely caused by the parallel slopes of 
the G-V curves, which indicate similar voltage sensitivi- 
ties of activation. 

In the presence of Ca 2+ , SRKR also decreases channel 
activity, at 1 and 100 uM, but not at 10 pM Ca 2+ . Thus, as 
well as the Ca 2+ -independent effect described above, 
SRKR also has Ca 2+ -dependent effects. Changes in sev- 
eral of the many parameters that define Ca 2+ - and volt- 
age-dependent gating could underlie these effects. 
Although the precise mechanism is not addressed in 
the experiments performed in this study, one potential 



explanation may be that BK SRK r gating is more sensitive 
to Ca 2+ than BKo gating. Supporting this idea, the de- 
crease in V1/2 with increasing Ca 2+ is greater for BK SRK r 
compared with BK). From zero Ca 2+ , the Vi/ 2 shifts are 
45, 188, and 241 mVfor BK SRK r, compared with 29, 142, 
and 211 mVfor BKo at 1, 10, and 100 pM Ca 2+ (Fig. 2 C 
and Table SI). Therefore, SRKR would be expected to 
inhibit gating via a Ca 2+ -independent mechanism but 
conversely increase the ability of Ca 2+ to gate the chan- 
nel. At 10 pM Ca 2+ , the lack of difference in the V1/2 of 
BKsrkr compared with BKo could be explained if the en- 
hanced Ca 2+ -dependent gating overrode the decreased 
Ca 2+ -independent gating of BK S rkr. This possibility is 
suggested by the observation that the range where no 
differences were observed in the V1/2 between BKo and 
BKsrkr (>1 and <100 pM) is near the Ca 2+ affinities esti- 
mated for both the RCK1 and the Ca 2+ -bowl sites (~4 pM; 
Bao etal., 2002). Thus, under conditions where the dif- 
ferences in the ability of Ca 2+ to gate the two variants are 
expected to be small (in the absence, at low, and at satu- 
rating Ca 2+ ) , the suppressive effect of SRKR would pre- 
vail, while the balance between the two gating processes 
would be altered at intermediate Ca 2+ (10 pM) . 

In addition to the complex effects resulting from the 
addition of the SRKR exon, we found that coexpression 
of the predominant brain BK channel accessory sub- 
unit, (34, had a range of effects on BKo and BK S rkr. As 
has been reported previously for other BK channel iso- 
forms, (34 slowed both activation and deactivation kinet- 
ics (Behrens et al., 2000; Brenner et al., 2000a). In 
addition, (34 altered BKo and BK S rkr G-V relationships 
in a Ca 2+ -dependent manner. Ca 2+ -dependent changes 
in (34 effects have been observed previously for other 
BK channel variants (Brenner et al., 2000a, Lippiat 
et al., 2003; Ha et al., 2004; Wang et al., 2006, 2009; 
Petrik and Brenner, 2007) . At low Ca 2+ , (34 inhibits BK 
currents, whereas at high Ca 2+ , it increases BK currents. 
A qualitatively similar process is occurring with BKo 
and BKsrkr, with a large rightward shift in G-V curves at 
10 pM Ca 2+ and small rightward (BK S rkr) or leftward 
(BKo) shifts at 100 pM Ca 2+ . Previous work has shown 
that (34 acts to decrease channel activation by increasing 
the energetic barrier for intrinsic pore opening and de- 
creasing coupling between the voltage sensors and Ca 2+ 
binding, and conversely, (34 acts to increase channel 
activation by increasing the coupling between Ca 2+ 
binding and pore opening, and stabilizing voltage-sen- 
sor activation (Wang et al., 2006; Contreras et al., 2012) . 
Given the broadly similar relationships between V1/2 
changes and Ca 2+ from previous work and our data, it is 
reasonable to think that (34 is impacting BKo and BK S rkr 
function via the same mechanisms, with the differences 
between the variants attributable to shifts in the Ca 2+ 
concentration at which (34 switches its net effect on the 
G-V relationship. This Ca 2+ concentration would be ex- 
pected to be higher for BK S rkr than BKq. However, the 
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relative interaction of (34's effects on Ca~ + dependence, 
activation, and deactivation kinetics remains to be de- 
termined for the observed differences between BKo and 
BK-srkr AP-induced currents. 
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